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NOTATIONS 


REPRODUCIBII.nY OF THE 
ORIGINAL PAGE IS POOR 


A - Area of a chamber section 

Aj^ _ - Area of a section occupied by burnt gases 

ApQ - Dome area 

Aj. - Area of a section defined by the flame front 

A. - Area of a jet section 

d 

A^ - Area wetted by flow per unit length 
Ajyj - Maximum chamber area 

A^ - Area occupied by flame at stabilization level 
Aqdo “ Area of kO’le.s, >p>>ac.ed in tkeddome 
Aqj^ - Area of cooling slots 

^OTF " air-injection ho'les.o.n the iflame rtuhe 

Arpp - Area of flame tube 

Ai ) 

A2 ^ - Constants involved in the combustion laws 

A3 ) 

B - Spalding transfer IjuncittEonn' 

Cjj - Coefficient of flow 

C-c - Coefficient of friction 

F 

G - Coefficient of drag 

C - Specific heat at constant pressure 

c 

d - Jet 'diameter - f - ac 

f - Characteristic parameter of impact pressure in the chamber 

P - Transfer parameter in two-phase flows 

I - Turbulence scale 

0 

m - Mass flow 

M - Mach number 

M - Molar mass 

P - Pressure 
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- Radius of curvature 
~ Mean radius of drops 

- Reynolds number relative to drops 

- Time 

“ Combustion time in turbulent flow 

- Dwell time 

- Temperature 

- Aerodynamic resistance 

- Rlow velocity 

~ Eddy ve-l‘0c'Lla;$"- ^ns 

- Laminar combustion rate* 

- Turbulent combustion velocity 

- Plow velOG'ity . \i 

- Drop velocity 

- Reduced length of chamber 

- Degree of progress of the reaction 

- Angle of jets 

- Temperature ratio before and after combustion 

- Specific heat ratio 

- Coefficient of thermal conductivity of the gas 

- Coefficient of viscosity of the gas 

- Density 

- Characteristic time of the vaporization of drops 


Subscripts 

A - Air 

b -BBurnt gases 

0 - Burning gas 

e - Intake 

1 - con dit ion s 

j - Jet, chamber radius mark 
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K - Kerosene 

g - Gas 

f - Fresh gas 

f - End of vaporization 

m - Mixture 

o - Initial conditions 

s - Flow point in a flow line 
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PUNDAMEHTAL ASPECTS OP THE AERODYNAMICS OF TURBOJET 

ENGINE COMBUSTORS 

M. Barr ere 

Scientific Director of Energetics at O.N.E..R.A- 


1. Introduction 


/55 


The improvement in turbojet combustion chambers has become 
increasingly difficult as a result of the restrictions imposed on 
the builder; the latter must therefore seek a compromise among 
the different trends [1]. The purpose of this report is to give 
an account of certain basic studies that should help the builder 
to achieve this compromise leading to the optimal solution for a 
given application. These studies are limited here to the aero- 
dynamics of the chamber'and affects both the main chamber and the 
af t erhurheni ch'amb er . 

Let us first examine the different qualities that modern 
combustion chambers must exhibit [2], [3]. 

- The first of these qualities is combustion efficiency, 
which must be as high as possible; excellent results have been 
obtained in the case of the main burner, but progress must still 
be made in .afterburner: 'chambers 

- The main combustion chamber must be as short as possible, 
while retaining the highest fuel efficiency. This minimum length 
is basic, since it reduces the .length of the shaft connecting the 
compressor and the turbine. For hre'a-S'crns of space, the length of 
the afterburner chamb-ers must also be a& s'Kort 'S&''p#sslble . 


* Numbers in the margin indicate pagination in the foreign text . 
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- It is equally important to reduce the pressure lasrsj in the 
chamber; it is certain that simplification in combustor design 

is a move in this direction ^ but one must not forget that this 
pressure 3i®’&S' is proportional to the square of the '■ 
intakeJMaCh number , A very careful aerodynamic analysis is neces- 
sary, however, if we wish to reduce pressure I’osjgies. 

- It is also necessary to design a cooling system for the 
walls that is as effective as possible in reducing the flow of 
air required for the protection thereof. This limitation is .'dif- 
ficult to comply with due to the increase in the richness of the 
chamber resulting from an increase in'_ compressionj’rati'o, inlet 
and' exhaust temperature of the main chamber. 

- In order to limit the temperature at the turbine inlet, it 
is necessary to dilute the burnt gases with colder air. The 
problems that thus arise are linked to homogeneity of temperature 

and flow compositionsat the turbine inlet. As a result of the 756 

intermittency in the structure of the flow, the turbine blades 
receive successively hot and cold bursts, but given the thermal 
inertia of the blades, the temperature of the material remains 
virtually constant; there is however a possibility of chemical 
reaction on the surface as a result of the variation in composi- 
tion that might contribute to a faster destruction of the blade. 

It is therefore necessary ,7lnsbfa~r ~^"po~ssi ble , t o p rovide a flqw\ 
ih^hbm^’erTeous as possible both in temperature and composition at 
the turbinef inOietbuibiiTc . 


- The builder is furthermore required to limit the emission 
of pollutants (solid particles and the production of nitrogen 
oxides). This emission is above all a function of temperature.., 
level and organization of combustion. 
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Though no attempt is made at a complete coverage of the sub- 
jects that serve as guidelines for the current research programs, 
we would like to underline in th course of this report the impor- 
tance of aerodynamics in the different problems involved in the 
improvement of a , combustion chamber. I 

The plan we have adopted is the following: 

-The first part is devoted to the main tombustdr listi.ng pa r-’ 
ticularly the aerodynamic problems relating to the preparation of 
the mixture in the recirculation zone where primary combustion 
takes place, and in the secondary combustion zone followed by the 
dilution zone in order to limit the turbine inlet temp'erajfure; 
the chamber ensemble is - studied to determine the .press.ure 
drop between inlet and exhaust; this pressure drop”^ depends not 
only on the velocity at the c,p.mbUstion chamber inlet but on the ' 
Qrganization„of the^wfldw as well. 

-The second part deals with the aerodynamics of the x *' -'-1 
chambers; a model providing for an evaluation of the development 
of pressure, temperature, velocity, and efficiency of combustion 
throughout the chamber has been developed. 

2 . Organization of Flow in the ’Ylain Chamber and Development of 
Combustion 


Figure 1 is a diagram of a main annular chamber of a turbojet. 
The development of combustion can be set forth diagrammatic ally 
in four stages : 

-Firstly, fuel Injection: in the example selected, the 

injection device, comprises an entire system whose role it is to 
introduce the fuel in the chamber in the form of the finest pos- 
sible drops, previously mixed with air. The injection device 
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used ISj for example, the S .N.E,C .M. A . device consisting of a 
T-shaped tube through which fuel drops and -air circulate at a 
high richness, of the order of 7 to 8, in order to avoid all com- 
bustion in th tube. 

-We then observe a recirculating flow zone, the place of 
combustion at a virtually stoichiometric ratio; the fuel is 
injected by the T-shaped tube into this vortex. 

-In order to complete combustion, the secondary combustion 
zone, which consists of the incompletely burned gases leaving the 
vortex region and of air injected through the wall, is situated 
downstream of the vortex. 

-Lastly, the final stage is the zone of dilution of burnt 
gases with air so as to bring the flow at the combustion chamber 
exhaust to a temperature consistent with turbine blading behavior. 

2.1. Injection Device 

\ 

The injection device plays a dual role; it must not only 

ensure a suitable supply of fuel, but it must, above all, provide 

- > 

for an atomized structure which, when mixed with the air, is the 
most likely to burn with proper combustion efficiency for a minimum 
chamber length, thus producing a minimum of pollutants, the ideal 
solution being reached with a gaseous fuel previously mixed with 
a portion of the air. 

In order to achieve this end, it is necessary to design a 
device thai; prevaporizes the fuel drops and produces a mixture 
with air while avoiding ignition of the medium. 

Pron the aerodynamic standpoint, it is therefore necessary to 
determine fast atomizing and vaporization conditions of the drops. 



The spray depends on the injection system used, and on flow 
conditions, especially on the difference in velocity between the 
jet and the air. In what follows, we assume the the structure 
of the atomized jet (mean diameter and distribution function of 
drops) is known. The study of the vaporization of fuel' drops in 
a flow is made on a theoretical plane from equations of conserva- 
tion of mass, of momentum and of energy established for the gas 
and the drops! 


The reduction of the mean radius r of the drop as a function 
of time is given by 'the relationship: 






v2__ 


where ro is the Initial radius of the drop and o = ^ 

this time constant 0' depends on the thermal conductivity of the 
gas of ■*^he density of the fuel 9 ^, of the specific heat of 

the gas Cp, and of the Spalding transfer function: 

n IC,jT7-:-T.) 1 
“ 1 

where T is the temperature of the gas surrounding the drop, T^ 
the vaporization temperature, AE^ the energy necessary to trans- 

from the fuel into gas. By introducing Prandtl number P, = , 

f a 


time constant 0 is written: 


0 ~ p r" — — - * ' 


The solution of the system formed by the equations of conser- /57 
tlon leads to the appearance of non-dimensional parameters that are 
representative of the development of the drops, and thus of the 
length of vaporization: 


-first, an overall transfer function 


9 Pr 9 

^'-2/r,Ti^]T)=2'p)rl 
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Fig. 1. Turbojet annular combustion chamber 


in which the Prandtl number and the Spalding number are involved; 
-then the Initial Reynolds number 
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Pf7 t O 1^0 

“ “ [iff 


Where Is the velocity of the gas and the velocity of the 
drops at the inlet plane; 


V 

-finally, the ratio (p = of drop velocity at the inlet to 
gas velocity at the end of vaporization. 


Length L of the combustion chamber Is represented by the non- /58 
dimensional parameter 

0M . 

J 

TEe pFactl!c¥i 'value's'~~¥f^F' ra nge from 0.5 to 1, and for 4) they 
range from 0.5 to 3. 



Pig. 2. Length of the vaporiza- 
tion tube. Effect of R 

eo 


One example of the results 
obtained is given in figure 2 
where is plotted as a func- 
tion of the Reynolds number 
for different values of (ji ‘and 
for P = 0.5; it is noted, 
particularly, that an increase 
in the Reynolds number reduces 
the length of the tube, while 
the 4» ratio is high. The 
effect of F about <j) = 1 
remains slight. This para- 
meter study provides indica- 
tions of the best definition 
of the vaporization device. 

It is specifically necessary 
to produce a homogeneous mix- 
ture of fuel and air, and to 
avoid' the vortices that run 
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the risk of becoming ignition zones. The problem, of ignition has 
been approached experimentally heretofore in order to determine 
the richness limit beyond which there is no ignition; an accelera- 
ted vaporization of the fuel results in a drop in temperature of 
the flow and thus avoids ignition. 

2.2. Recirculation Zone 

Side injection of air produces a turbulence, and this recir- 
culation zone is the site of a primary combustion with an air/fuel 
mixture in virtually stoichiometric proportions; the aerodymanlc 
study of this zone is complex, the important thing is the deter- 
mination of the air fraction entering the nucleus. In the pre- 
sent state of the art, this zone is likened to a Longwell homo- 
geneous reaction engine, whose op erat ion _i^' given by the balance of., 
species and the balance of energy. 

Let us take a chemical equation of the second order to char- 
acterize combustion of the form: 

K + A P 

Fuel + Air Combustion products 

Continuity of species P in the recirculation zone is given, 
in the steady state, by: ~ where is the burnt discharge 

leaving the turbulence, V the volume of the recirculation zone, 
and the mass rate of production of P due to chemical reaction. 
The preceding reaction leads to a value of of the form: 

where is the molar mass of P, and those of K and A, B the 
frequency factor, T the activation temperature, p and p the 

A ^ " 1 ^ 

density of the species A and K. Flow m leaving the recirculation 

• ^ P * * 

zone is equal to where mis the outgoing overall'' flow 

(burnt products and diluent) and p the mean density of the flow 
in the recirculation zone and at the outlet of this zone (perfectly 
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mixed Jet). By introducing gas dwall time t in the recirculation 


Pv 


m 


and defining the degree of advance of the 


zone defined hy t = 

s 

I 

reaction z = — -g—, the species P balance equation is written: 

B 


with stoichiometric mixture 








The equation of conservation of energy at steady speed is 
equal to: 

mCp (T — 'J'c) ’ 


Rise in temperature in the turbulent zone from to T is 

given by the input of energy by chemical reaction AH W , AH being 

C P P 

per unit mass of the energy released by a complete combustion; 
we assume that there are no losses. 


Assuming a maximum temperature released by combustion 

AH M 

— 2-R = _ T so that: 

Lae 

T — Tc - (T„ _ T,) /.R v'f (1 — 

or yet : _ _ 

T — Tc ^ (T.. — Tr) = 

we thus obtain a linear relationship between T and z. 

Let us introduce the production rate in the form: /^9 

C, ^ (1 _ zf ^ 2 ”' (1 — zf 

replaclng the exponent of 2 by a power of m degree of s, this 
approximation is valid for high values of the activation temper- 
ature (T^ > 12,000 K), being a characteristic time of the 
chemical reaction. 
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Pig. 3 - Effect of the imperfection 
of the mixture. 


On the same graph iFlg. 

33 iet'-us plot the varia- 
tlons of the I line 5^ 

and of the curve y = s'"(l— 
the solutions yielding the 
operation of the nucleus 
correspond to the inter- 
section of the straight 
line and of the curve 
(point A, for example). 

This simple theory enables 
us to explain the opera- 
tion of the recirculation 

zone. If 'dwell time t is 

s 

reduced, we arrive at 
point C that determines 
the blast conditions of 
flame; improvements have 


been made taking into account .mixing: time compared to dwell 

t 

time t . For a mixing time t = for example, we observe a 
s m pu 

deformation of the curve, the degree of progress of reaction 
is less, point A shifts to 'B for the same ratio of t -,/t . 


The aerodynamic study of this zone still has to be made, it 
should determine, specifically, the distribution of mass flow, 
momentum, and energy; certain attempts have been made by Spalding 
and Anasoulis using techniques of integration by finite differ- 
ences [5]. 


2.3. Secondary Uombustlon Zone 

Combustion efficiency in the recirculation zone remains low 
and ranges between 0.6 and 0 . 8 ; it is therefore necessary to 
complete combustion mixing the hot gases leaving the recirculation 
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zone with secondary air. This part of the combustor resembles a 
more or less stratified tubular reaction engine. I-n order to 
define this zone on an aerodynamic plane, we use the equation of 
one-dimensional balance in a manner similar to what would be done 
in the case of afterburner chambers, and this will be dealt with 
in greater detail. In a simple model, these two parts of the 
combustion chamber: recirculation zone and secondary combustion, 

zone may be likened to the juxtaposition of a homogeneous reaction 
engine and a tubular reaction engine . This coupling leads to a 
minimum dwell time, i.e., for a minimum combustor length if the 
homogeneous reaction engine is adjusted to the maximum rate of 
chemical production [2]. 


The flow diagram is generally the following: since the com- 
bustor is assumed to be annular, the flow tube situated between 

radii R. and R -_7 is considered i at a steady velocity, the balan- 
ces os mass, momentum and energy are written on the x abscissa: 


dm.] diiic ‘ , . 

-1- II 


dmt\ 

d . , > ^ dmt 

= mjQi 


■ w;) -h 


diHt 

dx 


lJ>JH 




j-iii 


{Cp]-l Ty_i Tj;) 


+ 


dnit 

lU 


j.j+1 




Where m. is the flow tube rate, m the external excess flow, 
the turbulent exchange rate between two neighboring flow tubes, 
p static pressure, T. arrest temperature and u axial velocity; 

Q is the rate of energy released by combustion [6], 


2.4. Dilution Zone 


The role of the dilution zone is to raise the temperature 
at the turbine inlet to an acceptable level for blade protection* 


11 





• ' "U ft *«-'o*. 

“ - -|,| .* ' Exhaust £ 

I H ot "gases ' 

'.Dxiution tJi 

» ' *. 


J • 

\ 


air 


L^_ 
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It is therefore necessary to 
iproduce in the shortest pos- 
sible distance a homogeneous 
flow both in temperature 
and composition. The assem- 
bly used for this purpose is 
shown in figure 4; dilution 
air circulates throughout the 
combustor i it is used to cool 
the wall providing a protec- 
tive film close to the sur- 
face; dilution air is injec- 
ted into the hot flow through large-sized holes. Prom the 
aerodynamic standpoint j, it is important to knowj on the one handj 
the mean line of the Jet entering the combustor and, on the other 
handj from this mean line, the diffusion of cold air in the hot 
gases, so as to be able to predict the degree of homogeneity 
of flow., in the exhaust plane of the combustor. 


Pig. 4. Organization of combus 
tion in the dilution zone. 
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2.4.1. Jet Penetration [7]“ 


The Jet forms an angle 3^ with the main stream; it is subjec- 
ted to forces that impose the curvature: according to the perpen- 

dicular n to the mean line, the centrifugal force applied to the 
mean line by an element having a length ds and a cross section Aj 
expressed as 

is equalized by the aerodynamic resistance T^. 

Along the tangent, the inertial force. 

P/«^/ j;- A/rfs 

is equalized by the longitudinal component of drag. 
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Fig. 5.' Diagram of the dynamic 
equilibrium of the jet. 


The equations of the 
equalization of forces and 
the equation of continuity 
yield the development of the 
mean line; the system is 
written: 


1 pti.2 Ct . 


■ nd 

pit* Ct (Vf \ . 

h[--cosjj5in- 

Vjd^ = tVorfj. 


di’j 

ds 


?- 



theory 

Pig. 6. Penetration compared 
according to the initial 
orientation of the jet, 

(mean line) — (v . /u = 8.32) 

3 ^ 


velocity and injeejjipn velocity to 


The .^notations used are ^ . 
derived from figure 5. 

The solution of this 
system leads to the deter- 
I mination of the mean line. 

I 

These results are compared 
with the experimental 
results for different values 
of as well as with 
Ziegler's calculation (fig- 
6), the agreement between 
theory and experiment is 
satisfactory . 

It is therefore possible 
in the case of a wide range 
of injection angles, flow 
have an idea of the mean line. 
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Poor 



Pig. 7a. Overall appearance of a dilution jet and 

wall film. I- 



Fig. 7b. Study of flow lines. 


2.4.2. Effect of Dilution 

The problem of the effect of dilution is more complex because 
it is difficult to describe^ experiments conducted in a tub of 
water have made it possible to show this complexity: thanks to a 
process of visualization, let us examine the structure of the jetj 
figure 7a shows the overall appearance of the dilution jet 
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underscoring the structure by bursts of the stream. A study of 
the flow lines Is made possible by a materialization thereof with 
a very localized Injection; the helical structure of the fluid 
thread surrounding the main jet (fig. 7b) Is seen. In order to 
show the Incidence of the dilution jet on the cooling of the 
walls, the visualization of the wall injection, downstream of the 
jet, shows (fig- Ja.) a return of the flow towards the dilution /6l 
jet. These visualization processes thus make it possible to 
describe the flow better 



Fig. 8. Study of dilution 

Key: (a) Axis of the hole 

(b) Axis of the exhaust 

section 

(c) Isotherms at exhaust of 

a dissymmetrical com- 
bustor 


A study of dilution 
by a distribution of tem- 
perature in the exhaust 
section of dilution has 
made it possible to deter- 
mine the structure of 
this zone. Figure 8 shows 
the isotherms at the ex- 
haust of a dissymetrical 
combustor, thus showing 
the efficiency of the 
cooling film (lower part 
of the figure [100 and 
300°'C isotherms] ) ‘ and the 
efficiency of dilution 
with the isotherm at i|00°C 
along the axis pf the 
two-dimensional combustor. 


2.5. Evolution of Pressure All Along the Combustor 


An important optimization relating to the_^erodynamics of the 
combustor deals with the evaluation of pressure loss. In effect, 
the performances . of the turbojet are affected by pressure loss 
in the combustor. These losses are due to three main causes: 
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the supply of energy in a subsonic stream, the losses due to fric- 
tion and separation, the losses due to jet mixing. If we desig- 
nate the impact pressure at the compressor exhaust as and 

the impact pressure at the turbine inlet as we define, gen- 

erally two criteria that characterize this loss in pressure: 

f Pi-'i.' — PtfT f _ — Pin 

' “ Pi^c ’ Jr 


the dynamic reference pressure being equal to 


<lr 


i 

A Ml. r j 

— ; ^ j 
^ M j 


1 rft^r 


where is the total rate of flow passing through the combustion 
chamber, density at the compressor exhaust, and Aj^ the max- 
imum annular cross section of the combustion chamber. q thus /62 

refers to a fictitious __ flow ~ where the stead y'sta t e of the fluid 
is the one corresponding to conditions at the end of compression 

p and T for static pressure and temperature, and au. = y nRT 
‘^sc sc ^ ^ ^ sc scsc 

for the speed of sound. We define a reference Mach number M 

r 

such as 

1 

I ’’ Am * I 

SO that : 


•• ro 

fjr = Psc Aif 




The study of an optimum combustion chamber leads to a maximuraa 
reduction of fi or fz . 


The strategy to be followed is thus the following: 


Since the conditions defining the upstream flow are given. 


X . e . 


^sc’ 


sc- 


M 


and sinc^/the'energyxinvolved in the 


combustor is defined by the arrest temperature ratioiiT ..pm/T , 

'I'SliX. *Z.^S c 

it is a matter of determining the geometry of the combustion 
chamber leading to the minimum value of fz {or fz) • This geometry 

J 

f . ■ , f 




V 



must also^ respond to the principles of combustion organization 
discussed in preceding paragraphs. 



Key: (a) ED Diffuser inlet; (b) Kerosene; (c) Dome (ADO 

area); (d) AOTP Combustion chamber air-supply hole 
Srei;; (e) ET Turbine inlet; (f) ATF mean flame-tube 
area; (g) SC Compressor exhaust; (h) AODO Area of 
holes in dome; (i) AM Maximum area of combustion 
chamber; (j) AOR Area of cooling slots; (k) Pri- 
mary air; (1) Secondary air; (m) Dilution air; 

(o) Axis. 


It is therefore necessary to anticipate the effect on f: of 

tke geometry of the diffuser (fig. 9) linking the compressor 
exhaust to the combustion chamber, the geometry of the nozzle 
directing the flow to the right of the flame tube, the cross 
section of the holes ^qj^q arranged on the dome, the straight por- 
tion of the dome the mean straight cross section of the flame 

tube Aijip and the geometry of the flame tube, the cross section of 
the holed situated on the wall of the flame tube ensuring 

the supply of air to the primary and secondary combustion zones. 


and to the dilution zone, and the cross sections of the cooling 

slots . 

UK 

The value of or fz is derived from the equations of con- 
servation of flow velocity of momentum and energy. We are not 
listing these classic equations, rather we shall list the hypo- 
theses one is lead to formulate, and the important terms involved 
in these relationships. 


There is not tinuch difficulty in writing the equations that~ 
define flow in the diffuser or the variable-geometry annular flame 
tube; it is enough to take into account a term of the type: 


f*s <IA 

J./ 


between abscissas aii and xz to evaluate the variation in cross 
section and of a friction term: 


J" 2 f ‘ P 

where C issthe coefficient of friction equal in the case of com- 
r 

bustion chambers to C = 0.0035 + 0.264 R and A the sur- 

r e m 

face per unit length wet by the fluid. Likewise, in the equation /63 
of conservation of energy, we can take into account the energy 
released by the combustion reaction assuMng a composition in 
equilibrium at the end of combustion. 


These difficulties in the evaluation of fi„ai^d fz are derived 
particularly from the holes arranged on the flame tube that dis- 
tribute air in the combustion chamber. Let us examine the effect 
of a hole. Flow in the annular space is divided into two parts, 
one relating to tracking the flow in the annular space, and the 
other to the _f low_ve3:o: c'i't'y./ through the holes: 


(a) Flow of the first part, which develops in the annular 
space, suffers an abrupt arrest due to the deviation of the flow 
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in^iithea hole, the loss of impact-rpr assure on either' side of the 
annular space is equal to: 



1 and 2 are the conditions in the annular flow before and after the 
hole, i.e., after deviation of the flow, mo is the deviated flow 

• ♦ f 

passing through the hole mi = m^o + mz 3 a and b are two constants, 
and a a function of the downstream Mach number. 

L = 4- 0,242 ' ; 

(b) rFlow m'o passing through tbe hol^^'can~bnly'be linked to the 

drop in pressure through the hole i^ffflow coefficient is known] 
this coefficient depends not only on the velocity through the hole, 
but also on the >pressure level in the flame tube.; value of{ for 

holes osfivarious shapes are currently available [8]-. 

(c) The jet enters the sflewm that occupies the flame tube; 
it is then necessary to know the angle of the jet with respecto to 
the mean flow in the flame tube, jet penetration and diffusion 
thereof, and its mixing with the main sfloTSr^i; as stated, the laws^ 
of penetration are starting to be well known, but /they are not so ^ 
,well known Inthe'h ase of the mixture; the results obtainedwwith fi 
and fz 3 assuming a perfect mixture or a stratified flow, are 
actually only slightly different, so that a quasi-instantanieous 
mixture can be assumed for this calculation without any major 
error . 

The minimum values of fmand fz are obtained by a parametric 
sweep. The main results are the following: 

-The diffused’- inlet Mach number and the reference Mach 
numberIsM , which is proportional thereto, have a substantial 
effect on the overall pressure drop in the combustion chamber. 
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as shown in figure 10 (the 
coefficient of proportionalty 
between and M being the 

ratio of ti^e maximum chamber 
area to the diffuser inlet 
area) . 

j Since the angle of the, j 

Idif fuser is d^ ermined by the 

separations it is not possible 

to reduce the M and thus 

<r 

increase the maximum area 

of the chamber without aepror-s- - 

hihii^iye lenghtening of the 

diffuser. M is therefore 
r 

fixed by the length of the 
diffuser selected. Different 
expedients are currently used 
in order to increase the angle, 
the solution that appears to be the most promising consists in 
sucking the boundary layer in the diffuser, the suctioned flow 
is of the' .order to 3 to H%; this suction also has the effect of 
standardizing the velocity. 

It is then a matter of selecting the gejametry of the flame 
tube that will lead to a minimum value of /a. This geometry can 
be determined by two coaxial cylinders cor by a divergent annu-lS-l’ 
tube. In each of these cases (fig. 11), we have an optimum value 
for the straight portion of dome this optiumumvvalueod'ep'en-- 

ding on the arrest temperature ratio between the combustion cham- 
ber exhaust and inlet. The lowest pressure Ql-oss is obtained with 
the divergent flame tube, fz Increases when the arrest temper- 
ature ratio increases . 
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Fig. 11. Pressure drop as a func- 
tion of the straight portion of the 
dome in the case of a flame tube 
c'dhsisting of two coaxial cylinders 
and of a divergent annular one [8]. 


this effect is generally slight; the 
reduces the value of fz considerable 

0.4 to 11 for A 


OTP^^M " 


The optimum value of 
Apo more generally, that 
of A^p the characterizes the 
area of a mean straight por- 
tion of the flame tube also 
depends upon the distribu- 
tion of the openings in the 
dome (AQpQ/^jy[) • increase 

of increases the 

value of App/Aj^ correspond- 
ing to the minimum of fi. 

In order to fix concepts, 
and for an arrest tempera- 
ture ratio of 2, is 

equal to 0.55 for A^j^^/AM 
= 0, A^p/Aj^ =0.6 for 

for 0.2. 


Likewise, the ratio 

that characterizes 

U 1 p M 

the flame openings ehangeS 
the optimum value of A^p/Aj^^; 
increase of A^^p/A^^, however, 

; fz going from l8 for = 

case M^= 0.63,i A^^^/Aj^j = 0.1, 


and - 2. 


The increase in the number of holes, for a constant value of 
AQipp/Aj^ reduces this criterion slightly. 


The cooling holes, characterized by their A^p area are' also 
involved in the calculation of fzi we note, generally, that fz 
decreases linearly when the ratio Aq^/Aq^PP increases. 


/64 
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The study of the overall aerodynamics of the combustion cham- 
ber is therefore very important; specifically, it permits, in the 
case of given inlet conditions, to determine the optimum geometry 
of the combustion chamber. 

3 . lAerodynamlcs oTf Afterburners 

This burnereis situated at the exhaust of the turbine, inlet 

flow is therefore rapid (at a ‘velocity o~f.the Brder of 100' m..sec ^ 

the static temperature is high (of the’ order of 1000°K) , and its 
composition is that of a combustion (kerosene-air) with a richness 
of the order of 0.2, i.e., with a substantial excess of air. The 
kerosene is injected upstream of the afterburner and, due to the 
temperature level encountered and the velocity, it is possi'fi'le td i 
assume thatttlh'e kerosene isi premixed with the gas emitted by the ' 
turblneet( homogeneous-type flame) downstream of the stabilization 
zone. Actually, there is a stratification of the* mixture ratio at 
the level of the injection system that modifies the composition'' of 
flow, but this stratification is neglected in the aerodynamic study. 

The flame is sta- 
bilized by annular and 
radial obstacles, gener- 
ally in the form of a 
trough as shown in fig- 
ure 12, Combustion 
takes place as follows: 


(a) Near the sta-' 
bilizers we observe a 
recirculation zone that 
acts as a homogeneous 
.- reaction engine, the' length I 
'of this' zone depends 1 


K 



in an afterburner. 


Key: (a) Stabilizer; (b) Flame zone; 

(c) End of combustion; (d) Rene 
circulation zone; (e) jet zone; 

(f) Interpenetration zone; 

( g) Homogeneous combustion zone; 
[ThJ Heterogeneous combustion. 
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on the velocity of the gas at the stabilizer level and on a chem- 
ical time that is a function of the mixture ratio. 

(b) The turbulent propagation of combustion, with stratified 
flow of fresh gases and burnt gases is located after this recir- 
culation zone, it is the Jet zone. 

Cc) Farther downstream, the fresh gas zone disappears and the 
flames that propagate after the stabilizers are mixed, it is the 
interpenetration zone . 

(d) Farther:’ downstream combustion continues in a homogeneous 
combustion zone; generally, and so as not to lengthen the combus- 
tion chamber, the homogeneous combustion zone is restricted to a 
value of an efficiency value of 0.8 i 0.9'SOthat the homogeneous 
combustion zone can be neglected. 

The important thing is to determine, for each section, combus- 
tion efficiency and the development of pressure, velocity, and 
temperature, so as to be in a position to calculat the heat trans- 
fer rate at the wall and the concentration of p^o.llut ants . "The 
model developed by O.N.E.R.A. [9] neglects the circulation zone, 
the most important part of the combustion zone consisting of the 
turbulent propagation of the flame as shown in figure 13 . 



O A 


Xa 

Fig. 13 . Diagram of the afterburner 
Key: (a) Fresh gas zone; (b) Combustion zone; (c) Burnt gas zone. 
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Since , flow veloc ity_ — high with nespect to the turbulent 
combustion rate, the combustion zone (index c) is very long and 
thick and it separates the fresh gas zone (index f) from the burnt 
gas zone (index b). The way the system of equations governing 
flow is written depends on the development of combustion; three 
zones are marked: 


-The zone downstream from the -portion where *the burnt gases 
^re situated '( zone (^ ) , I 

-The zone downstream from the portion where the fresh gas 
ends (zone , 

-and, between these two, zone . 


Let us assume that for an abscissa of a: = 0, static pressure /65 
is ~p^', "ficfw veTocXty , static temperature T^, density P^, and 

the area occupied by the flame (recirculation zone) A^, the overall 
continuity equation for zone 2, for example for the abscissa 


X = X „ J 


Po^o (A — A(j) — p/**/ (A — A/) -}- pftWsAs -j- 


fA/ 

+ pcUcrfA 

J At, 


A is -the cross section of airflow ^ velocity is assumed to be 
uniform in the fresh gas zone and in the buiiit gas zone. 


Likewise, the equation representing momentum becomes : 


iJoA -f Potig (A — Afl) = pA 4- (A — A/) 

I + PbitlAt + J ^ 

To these two equations we must add the equation of conserva- 
tion of energy for fresh gases defining T^ and u^: 

1 CpTo + 2^ = CpT/ 

and the equation of conservation of energy for burnt gases: 
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where hh is the energy brought into play by combustion. 


The gas Is assumed to be perfect ^ so that: 

p p_T 

po PoTo 

and the development of p and T Is assumed to be isentropic: 



In order to obtain the aerodynamic fields it is first neces- 
sary to characterize turbulent combustion; two magnitudes of tur- 
bulent combustion phenomena are used, the first is the turbulent 
combustion rate which will define the geometry of the area 

separating the fresh gas zone from the combustion zone and, par- 
ticularly the inclination of the flame front with respect to the 
mean flow is given by the quotient of axial flow velocity u by the 
turbulent combustion rate the second characteristic magni- 

tude of the combustion phenomena is the time necessary for com- 
bustion, which will determine the thickness of the flame. If, 
for example, we follow a flow line F, E, G, K, the end of com- 
bustion, l.e.,,point G, is defined by , s being the 

' ' J P K 

distance travelled along a flow line on which, by virtue of the 
slight inclination of EG with respect to the axis, the local flow 
yalocity resembles axial velocity: u.. It may be assumed that 
turbul nt state combustion occurs Inbursts whose initial dimen- 
sion depends on the scale of the turbulence the combustion 

rate of these bursts u-^ is equal to 

uji = itfci. "i » 


where is the laminar combustion rate, I the dimension of the 
burst at moment t and is the eddy velocity, velocity 




Fig, l4. Simplified diagram of'the combustion chamber 


Key: (a) Confluence secti on; (b) Flame fron t; (c) Fresh mixture; 

(d) Combustion; (e) gremixed alr/metha ne; (f) Flow line; 

(g) Burnt gases; (h) Internal hot gases. 



= — dl/dt so that after integration combustion time is 
written: 



Aj is a constant (Ai = =; 4.3) this relationship has been duly 
verified by experiment s, combustion time depends on three para- 
meters, laminar flow velocity u-,ri thals'ca'le x©':f'Dtiunbul-enc,e cX, 

A ' foL ’ o o 

and eddy^ve'liOclt'yx,. When turbulence is substantial (w' ^ 

this is the case of an afterburner chamber, the turbulent com- 
bustion rate is given by 
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Fig. 16. Distribution of pressure along'^fcbeaa2;3;s . 


Key: (a) Pressure (abs. bar); (b) Experiment; (c) Theory. 



Pig. 17 . .Distribution of temperature. 

Key: (a) Flame front; (b) End of combustion zone; (c) theoretical 

curves; (d) Experimental curves. 
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i UfrT — AjPiibi, 



where 6 Is the temperature-combustion ratio 3 kz and A3 are con- 
stants, when combustion is confined Az = 1 and A 3 = 2; this last 
magnitude A3 = 2 stresses the role of combustion as a generator ‘ /67 
of turbulence - 


This model was first verified by means of a simplified 
experiment, the diagram of which is shown in figure l 4 ; the air- 
methane premix is initiated by a parallel flow of hot gases, the 
turbulent flame develops from this point. The velocity profile 
is the one in figure 15; a linear distribution of velocity is 
assumed in the flame zone. 


The only parameter that is not known is the e'ddy veio- - 
city u which has to be introduced in the expression of 
turbulent combustion velocity and combustion time. The comparison 
of theory and experiment deals primarily with the development of 
pressure throughout the combustion chamber; by adjusting the 
value of M ' we have, as shown in figure 16, a good agreement 
between theory and experiment. The geometry of the flame front 
is quite well described by this model, and so is the development 
of the isotherms (.fig. 17 ). This model has been applied to 
actual revolution combustion chamber, and the results obtained 
also agree ^Jiith experiments, verification deals wit-h the meanag-^ 
combustion efficiency. 


We thus have means for describing the development of com- 
bustion in an afterburner chamber (pressure, temperature, velocity); 
this model is based on the knowledge of the turbulent dddy^veioen 
city tlitl also assumes that the mixture is Initially homogeneous; 


bv means of 


Htn'no* of t hp nrcp-re m- 


28 





it is possible to introduce a certain distribution of the rich- 
ness downstream. 

4. Conclusion 


Knowledge of the structure of flow in the combustion chamber 
is of great importance if we wish to advance in this field. In 
closing^ we shall stress certain specific points. 

(1) The fuel-air mixture zone is essential in the operation 

r » 

of a short combustion chamber or a chamber that is only slightly 
polluting. The analysis of this zone is based on the study of 
two-phase flows, and particularly on the study of the vapor Iza- 
tioncof the drops in a hot turbulent flow, which is an area for 
which we have only few practical results as yet. 

(2) Combustion starts in the recirculation zone. The aero- 
dynamics of this zone -that is likened 'to a homogeneous reaction I 
engine a^ very poorly known and basic research in this area is 
necessary. 

(3) The Idllution zone li mits the temperature downstream 
from the combustor ; 'even- there , while it is possible to duly 
represent the development of the mean flow, there are few results 
insofar as the turbulent diffusion of Jets and the degree of 
homogeneity at the ^chamber exhaust are cohcernej d. 

(4) In order to produce a combustion chamber with a minimum 
pressure loss, it is necessary to prepare a realistic model of 
this ccombustion chamber. This modeliis, for the time being, 

f 

simplistic, and a great effort should be made to Improve the 
situation 
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(5) The aerodynamics of the afterburner chamber are based 
on a better knowledge of the development of turbulent combustion; 
In order to advance ^ it is therefore urgent ailso thavievmor.erdaiSata 
on turbulent combustion; it is particularly necessary to have a 
better knowledge of the levels and scales of turbulence upstream 
and in the flame . 
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